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We develop an approach to calculate the admittance of effectively one-dimensional open quantum 
systems in random phase approximation. The stationary, unperturbed system is described within 
the Landauer-Biittiker formalism taking into account the Coulomb interaction in the Hartree ap- 
proximation. The dynamic changes in the effective potential are calculated microscopically from 
the charge-charge correlation function resulting from the stationary scattering states. We provide 

■ explicit RPA-expressions for the quantum admittance. As a first example the case of a quantum ca- 

pacitor is considered where we can derive a small-frequency expansion for the admittance which lends 

, itself to an experimental testing of the theory. A comparison of the low-frequency expansion with 

the complete RPA-expression shows that for a quantum capacitor a simple classical equivalent cir- 
cuit with frequency-independent elements does not describe satisfactorily the quantum-admittance 
with increasing the frequency. 
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; I. INTRODUCTION 

The ac-transport properties of nearly ballistic devices are interesting from the point of view of basic research as 
well as from the point of view of technological applications. In basic research ac-transport can provide valuable extra 
information to stationary transport. In applications impedances are of crucial importance for the lay-out of ac-circuits. 
In the past years a number of methods for the description of time-dependent transport phenomena in mesoscopic 

\ systems have been proposed (for a recent review see Ref. ) . Among these approaches are techniques based on non- 
equilibrium Greens functions^Si^*^ Wigner functions)^ bosonization schemesfi and phenomenological considerationsi^ 
If there arc only small ac-fields it seems promising to calculate the linear ac-response to a perturbation of the 
stationary system described in t he successfu l Landauer-Biittiker formalism (LBF) jgiifliiiiiSiiiii Such an approach has 

• been developed in Refs. PIT^Il7liall<i2(T| and it has been used in a number of applicationsiSiiSSiSiS^iS^iSSiSi Here 

', the response to an external potential is derived 'which prescribes the potentials ([/„) in the reservoirs only' (see Ref. 
|l5l|). The reservoirs carry the charge Qa so that the ac-potential leads to a perturbation in the Hamiltonian as given 

[ by Hi = X]q UaQa- This approach avoids the calculation of the time-dependent potential within the scattering area, 
i. e. outside the contact reservoirs. However, the knowledge of the microscopic potential in the scattering area is 
necessary to derive a formal response theory, beyond the often invoked spatially uniform electric field perturbation^ 
As pointed out in Refs. '27"29''30(| the appropriate response formalis m for the interacting electron system is the 
random phase approximation (RPA). In our previous papers Refs. I3J32^ we demonstrated the application of 
the complete RPA-schemc to open stationary systems described in the Landauer-Biittiker formalism. A complete 
RPA-scheme2^ requires as a central clement the calculation of the irreducible polarization Ilo{r,r' ,uj) from the self 

[ consistent scattering functions of the stationary system. As a second necessary ingredient for the implementation of 
the complete RPA-scheme, it was shown in Refs. [sills^ how the dynamic total potential in the scattering area can 
be determined microscopically using the calculated irreducible polarization and the Greens function for the Poisson 
equation with Dirichlct boundary conditions. 

After a formal derivation of our theoretical approach we derive in this paper explicit RPA-expressions for the 
frequency dependent impedance in a general two-terminal device under large dc bias. These expressions are evaluated 

' for the case of a quantum capacitor. In the limit w ^ an expansion of the admittance follows as given by 

Y^-iij{Yi+iu}Y2 + ...), (1) 

with real constants Yi and Y2 that can be calculated directly from the scattering functions of the stationary system. In 
numerical computations we determine the admittance of a MIS-type heterostructureS on which measurements of the 
static capacitance have already been made^ In order to propose an experimental test of our RPA-approach we first 
compute the drain-source- voltage dependence of the coefficients Yi and Y2 in the limit a; — > 0. For higher frequencies 
we find numerically, first, that the expansion in Eq. becomes invalid very quickly and, second, that an equivalent 
circuit with frequency independent R and C elements does not reflect correctly the ac-propcrties of the considered 
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system. Instead, there are pronounced and systematic deviations from the equivalent circuit behavior which should 
be testable in experiments as well. 



II. GENERAL ADMITTANCE FORMULA 
A. Stationary system 
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FIG. 1: Sketch of the potential energy along the growth direction in a planar two-terminal system. It can be seen that 
—eUsD = V2 + p.2 — Vi — fii, where jl2 and p.i denote the diflerence between the chemical potential and the bottom of the 
conduction band in the bulk contact. 



We consider a planar two-terminal system with a source contact {s = 1, z < zs) and a drain contact (s = 2, z > z^) 
under an external bias Usd ss shown schematically in Fig. ^ The material of the contacts can either be a strongly 
n-doped semiconductor or a metalization. In the mean field theory the effective potential energy is independent of the 
perpendicular coordinate r± = {x,y), V — V{z). Because of the effective screening the potential energy is constant 
in the bulk of the source contact, V{z < —d) = Vi, and in the bulk of the drain contact, V{z > d) ~ ¥2- We write for 
the wave functionS^ 



ip{r) 



^(,)^^^^^exp(zfe^rj 



(2) 



where k± = {71^271/ L^, ny2Tr / Ly) , and Uy are integer numbers, and A = L^Ly is the cell area of periodic boundary 
conditions in the perpendicular directions. It follows from the time-independent Schrodinger equation {Hq^ E)ip = 0, 
with Ho = -?iV(2m*)A + V{z), that 



d^ 

2m* dz"^ 



V{z) 



(3) 



with 



E = e 



2m* 



(4) 



The wave function in Eq. |(2Jl is defined by the set of quantum numbers (s, e, kj_) without considering the spin quantum 
number. Consistent with the Landauer-Biittiker formalism the ip^^\e,z) are scattering functions, where the index 
s defines the direction of incidence: the scattering functions incident from the source contact (s = 1) exhibit the 
asymptotics 



/27r 



y^-^ ( exp [iki{z + d)] + Sii{e) exp [—iki{z + d)], 
S'2i(e) e^p[ik2{z - d)], 



for z < —d 
for z > d. 



(5) 



where and S21 are elements of the energy dependent 2x2 scattering matrix^ 

1 2m* 



(6) 
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and d{e — Vi) is the step function. Analogous expressions hold for the drain-incident scattering functions ifj'^'^^e, z). 

As usual, the effective (total) potential energy contains an external part, T4a;t, and a second part, V^i, coming from 
the Coulomb interaction between the electrons. The external potential arises typically from different band offsets in 
the used materials or from fixed external charges like fully ionized impurities. The Coulomb interaction between the 
electrons is taken in the Hartree approximation so that Vei obeys the Poisson equation having as sources only the 
electronic charge density 

^V,iiz)^~-Piz), (7) 

CIZ Kig 

where Hs is the dielectric constant of the host material. In Appendix IXI we reproduce that the electron density 

P(^) = 2-^^/ d6 5.(6) V-^^) (6, z) In {l + exp(/3(A..- e))} (8) 

in a formal quantum statistical approach needed to formulate the linear response theory which we will describe in the 
next section. In Eq. © gs(e) = m* /[h^ks{e)] is the one-dimensional density of states, (3 = l/{kBT), and we included 
a factor of two to account for the spin degeneracy. 

B. Harmonic perturbation 

1. Random phase approximation 

We consider our system with an additional small ac bias SU superimposed to the source-drain bias so that 

UsD{t) = UsD+SUe-'^^+''^^\ (9) 

where 77 ^ 0, ry > 0, is an adiabatic turning-on parameter. Because of the good screening in the contacts the applied ac 
bias is assumed to lead to a dynamic potential perturbation S(f>{r, t) = S(j){r) cxp [—i{uj + iri)t\ fulfilling the boundary 
condition 5(j){x,y^z < —d) = and 6(f>{x,y,z > d) = SU. Then, the perturbation in the potential energy felt by 
electrons is 5V{r, t) = —e5(l>{r) exp(— 1(0; -t- iri)t) and the time dependent Hamiltonian becomes 

H = Hn+ J d\ p{r)SV{r, t), (10) 

where /5(r) is the particle density operator. Of course, Eq. (|10|l is only valid if retardation effects in the quantum system 
can be neglected, i. e. w ^ c/L, where L is the typical length of the device. For a quantum device with i = 10 nm we 
find c/L ~ lO^^Hz, which is on the upper limit of the C/y-radiation. The time-dependent perturbation in the electron 
density (induced density) Sp{r, t) is calculated outside the contacts in random phase approximation Assuming 
for the planar structure SV{r,t) = SV{z,t) we show in Appendix IbI that Sp{r,t) ~ Sp{z,t) ~ Sp{z) ex-p{—iujt), with 

Sp{\z\<d)= [ dz'no{z,z\u})SV{z'). (11) 
Here Ilo{z,z',uj) = EL'=i no''^ \z,z',uj), 



t M r°° F^^^'Mz z' e f'^ 

Iit\z,z',u) = \in. de d6'i-^-^lf^, 
Jv, Jv, e-e' + h{uj + 11]) 



[12) 



and 

1 + exp (/3(/is - e)) 



F(-')(z,z',6,e') = 2-^.g,(e)g,,(e')ln 



27Tph^ w;.. V y [ 1 -f exp {P{p,, - e')) j 
(^(^)(e,z))V(^')(e',z)(v(^')(e',z'))V<^^(e,;2')- (13) 
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Setting the integration limits in Eq. Hll(l we assume that there is no phase coherence of the wave functions between 

the contacts and the scattering area, IIq'"* '(z, \z'\ > d,uj) = 0. The perturbation of the effective potential outside the 
contacts which enters Eqs. pO|l and Hll|l is determined by the Poisson equation 

ASV{z) = Sp{z). (14) 

Kg 

In Eq. (|14|l we assume that there is no mobile charge in the interval —d<z<d other than that of the tunneling 
electrons. The solution of Eq. (|14|l obeying the boundary conditions dV{z < — c?) = and dV{z > d) = —eSU can be 
written as 

5V{z)=5Vo{z)+ f dz'vo{z,z')Sp{z'), (15) 
J~d 

with the homogeneous solution 

5Voi\z\ <d) = -eSU^. (16) 

In Eq. (|15|1 the symmetrical Green's function vo{z, z') = —{e^ /2ks)[\z — z'\ + zz' jd — d\ for the Poisson equation obeys 
the boundary condition^ vq{z = id, z') = 0. Using Eqs. 111|) and (|15|l one obtains an integral equation for the total 
potential as given by 

/d pd 
dz' dz"vo{z,z')IlQ{z',z\uo)5V{z"). (17) 
-d J-d 

We write down the inverse of this equation in a convenient discrctized form 

SV = {l-voUoy^SVo, (18) 

with z ^ Zi ~ —d + {i — l)Az, i ~ I...N + I, Az ~ 2d/N, and ^ oo so that J^^ dz Az X^ila^- Furthermore 
we define the iV + 1 x + 1-matrices Vq and IIo{llj) with (t'o)ij = Azvo{zi, Zj) and (n-o)ij — AzIlQ^Zi, Zj,uj) as well 
as the vectors SVq and 6V with {6Vo)i = 6Vo{zi) and {SV)i ~ SV{zi). The continuum limit of H18|l can be regained 

using the von Neumann theorem, (l — VqUo) = S^o(''^o-'^o)", and rewriting the obtained sums as integrals. We 
solve Eq. (|18|l numerically. Then after defining the vectors dp with {Sp)i = 5p(zi) and an analogous vector 5jz for 
the z-component of the particle current density one obtains from discretization of Eq. (|11|) 

5p = Tl^5V = Il^{l-v^Tl^y^5VQ, (19) 

and from Eq. ljB5|l 

6j,=tloSV = tl^{l-voUoY^5Vo, (20) 

where Hq is the current-density response function defined and evaluated in Appendix^ and Ilg is the corresponding 
matrix obtained after discretization. 

2. AC- admittance 

From the continuity equation {d/dz)6jz{z,t) = —{d/dt)5p{z,t) one obtains the relation 

5jz{z) = 5jz{-d) + f dz'SUz') = 6jz{~d) + iLuSQiz) = -51/ Ae + iu5Q{z), (21) 

J -d 

with 5j'^{z) ~ {d/dz)Sjz and SQ{z) = J^^dz'Sp{z'). From the boundary conditions for SV{z > d) = —eSU and 
6V{z < — d) — and the continuity of {d/dz)dV{z) it follows that the total induced charge 6Q = 5Q{d) vanishes. One 
then obtains Sjz{—d) = Sjz{d) = —SI/Ae, where SI is the induced electrical current provided by an external source 
and flowing through the device. The minus sign in the last step of Eq. (|21|l results from the current convention. The 
complex admittance is defined as usual hy Y = 6I/SU = —Ae5jz{—d)/6U. Applying Eqs. (|21|1 . H20|l . and (|16() one 
obtains an explicit expression for the admittance 

Y^^=Ae^W[U„{l-VQnoy'Wo, (22) 
where (Wo)i — Zi/2d + 0.5, and Wi is a unit vector {Wi)i — Su. 



III. QUANTUM CAPACITOR IN RANDOM PHASE APPROXIMATION 

A. General results 

We define a quantum capacitor through two conditions: First, there is basically no dc currents traversing the 
structure, 

V'(^'(e,z > d) === 0, V*^^(e,z < -d) = 0. (23) 
Second, the overlap of the right-incident and the left-incident scattering functions can be neglected, 

^W(e,z)^(2)(e',z)^0, (24) 

for all e,e' and z S We then find from Eqs. and ^ n|,^^^ = H^^^^ n|,^^^ = n[,^^^ = 0, and the 

induced current can be split into two independent parts, ~ ^J'z'' + with 

<5jt' {z) = dz'nf ) (z, z\ u:)W{z'). (25) 

J -d 

is) 

Since bji results exclusively from the source-incident scattering states for s = 1 or exclusively from the drain-incident 
scattering states for s = 2 we can write for each component a separate continuity equation, 

-*c.<5p(^)(z) + ^5jt^(z) = 0, (26) 

with 

Sp^'\z)= dz'Il'^''\z,z',u;)dViz'). (27) 
Integrating Eq. H26(l one obtains using Eq. H23() under another form, i.e. Sji^\d) = Sji'^\—d) = 0, 

i^- = -5j(i)(_d) = ic^ / dz [ dz'U^o^^\z,z',u;)6V{z'). 

6^ J-d J-d 

With the definition of the admittance and Eqs. p8|) and the equation H22|l reduces after discretization to 

Y = -e2Ait^AzWjn[|"'(l - voUay^Wo, (29) 

where Ho is eliminated. In Eq. H29l we define the + 1 x iV + 1-matrix (n[,''''),y = AzU''„"'\z,, Zj , uj) and the vector 
(W2). = 1. 

In Appendix [ni it is shown that for small frequencies an expansion 

4^^) {z, z', cj) = P^-'^ [z, z') + iujP['^ (z, z') (30) 

can be derived with real functions P^^\z,z') and p[^\z,z'). Inserting this expansion in Eq. ()29|l one obtains a 
low-frequency expansion for the admittance of a quantum capacitor as 

r w -w(ri + wFa). (31) 

Here the leading order coefficient 

Fi = e^AAzWa^Pfl'^ (l - voP^Y^W^, (32) 



(28) 



and the first correction 



Y2 = e^AAzW'^ 



p[^\l-voPa) ^ + P^^\l - v^P^) 't;oPilWo, (33) 



are real, where Pq = Po^'^ + Po^'\ Pi = Pf ^ + Pf \ {P^o%j = AzP^'\z,, z, ,u;), and (P^'^),, = AzP^'^z,, z,,u;) . 
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B. Numerical results for MIS-type nanostructure 



As a test structure for our theory we take a planar MIS (metal insulator semiconductor )-type GaAs/A^Gai-^^As 
heterostructure with a near back gate. This structure has been analyzed in experiments)^ the stationary system 
has been described theor etically within Hartree approximation^ and first calculations for the dynamic behavior are 
presented in Refs. |3ll32l |. We now calculate the quantum admittance Y{u)) according to Eq. (|29|l up to frequencies 
of lOOGHz. Using Eq. one can extract from the numerically (or experimentally) given data the parameters 
Yi = — lim^^^o Im[Y{LLi)]/LLi and Y2 = lim^^^o Re[Y {oj)] / ui'^ . It is then possible to recast Eq. H31(l in a normalized form 



Y 



'iu!{l + iuj), 



(34) 



with Lo = Lolu}^), LUo = Y1/Y2 and Y = YY2/Y1. This normalization allows us to collapse the calculated quantum 
admittance at all considered source-drain voltages into one graph which is presented in Fig. |2| It is immediately 
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FIG. 2: Quantum admittance for MIS-type nanostructure, in a normalized representation, for different static biases Usd, 
parameters corresponding to Fig. 1 of Ref. lioll . Dashed line represents Ysg{ui) for a classical RC-circuit with frequency- 
independent elements. The symbols represent the critical frequencies— up to which the approximation of the quantum result 
(i.e. Eq. (^J) with Y^g can be considered satisfactory. 

seen that the expansion in Eq. (|34|l only holds for (D ^ 0. For finite frequencies there are significant deviations from 
the value 1 for Re{Y)/Cp and from the value —1 for Im{Y)/Cj. To discuss these deviations we compare with the 
admittance Y^g of a classical equivalent circuit consisting of a frequency-independent resistance R and a frequency- 



independent capacitor C in scries 



17.41.42 



An inspection of the admittance Ysg{uj 0) of this circuit yields Yi = C 



and Y2 = RC . These formulae can be regarded as quantum mechanical expressions for the elements of the equivalent 
circuit. In the normalized form one then obtains 



(35) 



It is seen from Fig. |5]that Ysg{Lu) generally fails to describe the numerical admittance. The numerical results show 
for increasing source-drain voltages a systematic enhancement of the decrease in Re{Y) and the increase in Im{Y) as 
the frequency is increased. This finding docs not result in the classical equivalent circuit. 

In Fig. O (a) we represent the dependence of the coefficients Yi and Y2 on the working point bias Usd- The 
coefficient Yi is essentially identical with the low frequency limit of the dynamic capacitance plotted in Fig. 3 of 
Ref. Is^l- As demonstrated in Ref. [s^l this step in Yi is in good agreement with a step in the experimental 
capacitance curve which is caused by the formation of a two-dimensional electron gas within the quantum capacitor. 
The coefficient Y2 shows a general increase with increasing bias. As a characteristic feature it is seen that the step in 
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FIG. 3: a) Yi and Y2 as function of the working point Usd- b) The scahng frequency wo vs. Usd- 

the capacitance of Yi is accompanied by a small hump in ¥2- An inspection of the scaling frequency loq (Fig. El (b)) 
reveals a corresponding downward hump. 



IV. CONCLUSIONS 



We present a quantum mechanical model to calculate the admittance of effectively one-dimensional open quantum 
systems in random phase approximation. Explicit RPA-expressions for the quantum admittance of a general two- 
terminal system are derived. In the case of a quantum capacitor a small-frequency expansion can be obtained which 
lends itself to an experimental testing of the theory. A comparison of the low-frequency expansion with the complete 
RPA-expression shows that for a quantum capacitor a simple classical equivalent circuit with frequency-independent 
elements does not describe satisfactorily the quantum-admittance with increasing the frequency. 
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APPENDIX A: STATISTICAL OPERATOR FOR THE STATIONARY SYSTEM WITH FINITE 

SOURCE-DRAIN BIAS 

1. Scattering states as a complete orthonormal single particle basis 

As shown in Ref. [s^ the scattering states ^'''^^{e,z) constitute a complete orthonormal system, i. e. 

^ r&g,(e)(^(^)(e,z))*V^^)(e,2')=<5(^-^'), (Al) 

s=l,2 

and 

r dz U^'^^ie, z))* V^^'He', 2) = ^(e " V,)Sss'Sie - e')/gs{e). (A2) 

J —00 
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To eliminate the weight function (7s (e) we substitute for a given s in Eq. ijHJ e = h^k1/{2m*) + Vs = £s{ks) (see Eq. 
©). One obtains 

r dh{i;^'\eM,z)y yj^'^{es{h),z')=S{z-z'). (A3) 

1 O 



According to Eq. © we write in Eq. HA2p e' = {k'^,)'^ / {2m*)+Vs' = es'(fc^,) and using the identity (5ss"5(e—e')/gs(e) = 
Sss'^iks — fcs) it results that 

dzU^'\es{k,),z)y i,^''\e,,{k',,),z)^0{k,)S,,,6{k,~k'^). (A4) 

Since later we want to work in a number representation we introduce a /c-space discretization kg — > kj = jAk, 
j = 0,1,2...., so that S{ks — k'^) 6jj'/Ak. Furthermore, we define ^sj{z) — V Ak ^p'^'^^ (^es{kj), z) . With this 
definition an explicit asymptotic form of the source-incident scattering wave functions (s — 1) follows from Eq. ((SJ as 
given by 

I ( exp (ikj{z + d)) + S'ii(ei(fcj)) exp (^—ikj(z + d)), for z < ~d 

i^lA^)^^]^0{e,{k,)-V^))\ (A5) 

^ [ 52i(ei(fcj)) exp [i^k-] + ^{Vi - V2){z - for z>d. 

For the drain-incident scattering functions (s = 2) we find the asymptotics 



'Afc. 



,(e2(fc,)) exp (^-^^^ + ^J^{V2-V,){z + d)), 



for z < —d 



exp (— ifcj(2; — d)) + S'22(e2(fcj)) exp (ifcj(z — d)), for z > d. 
After discretization we thus write for ljA3|l 

2 oo 

5^5^^:,-(z)V'.,(/)='5(^-^'), (A7) 

s=l j=0 

and for HA4|I 

/•oo 

dz'ip*j{z)ips'j'iz) = Sss'Sjj'. (A8) 

In addition we introduce a complete orthonormal basis system 0i/(rj^) for the square- integrable function in /i^, 
/ dr±(f>l{r±)(l)^i{r±) = 5vu', and ^y4'Z{''''i-)4'y{i'' i.) = ^('"J. — '"'_l)- They are usually solutions of the time- 
independent Schrodinger equation in the lateral directions, [—h^/{2m*)A± + V±{r±) — i?^]0y(rj^) = 0. Then a 
complete orthonormal basis for the Hilbert space of the single particle quantum states is given by 

Va{r) r\ip^ >^ i'sjiz)Mr±), (A9) 

where a is the index-triple {sjv). We find the usual completeness relation 'Yl,a.f*a{''')'f = ^{i' ^ '''') and the 
orthonormality relation is given by J d^rip'^{r)ipa>{r) ~ ^aa'- 



2. Definition of the statistical operator 



Because the scattering states |(/3q, > constitute a discrete and complete orthonormal basis it is possible to introduce 
the creation- and annihilation operator cj, and Cq, respectively, with the usual anticommutation relations {c^ , c^/} = 
<5a,a', {cct, Cq'} = 0, and {cj,, c]j;} = 0. Based on the anticommutation relations one can formulate a particle number 
representation which will be described in the following. 

According to standard theory^ the field operators are given by 4'(r) — J2a'fia{''')ca and ^''{r) = X]q ¥'a('")^Q- 
The many-particle Hamiltonian of the stationary electron system can be written as 



Hn = 



d^r 



-—A + V^{r^) + Viz) 



*(r). 



(AlO) 
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where Ea = E'^ + €s{ki) is the energy of the single particle state \(pa >■ Here we use the notation Hq since in Sect. 
IllBI we consider a time dependent perturbation of this Hamiltonian. We now use as a statistical operator for the 
stationary system as given by 

po = ^ cxph/3(i/o - Mi^i - M2A^2)], (All) 

^0 

with the chemical potentials fix and fj,2 of source- and drain contact, respectively, the particle number operators 
= E^u ci,^Cs^^. with N = Ni + N2, and 

Za = Tr{cxp[-/3(Ho - MiiVi - ^i2N2]}■ (A12) 

The trace is done over all states of a Fock-space basis, constructed with the help of the single particle scattering states 
\ipa >■ The trace is easy to write in the occupation number representation 

00 

Zo = '^^<N;...nc,...\exp[~P^{Ea,-fi,,)cl,Ca']\N;...na---> 

N=0{na} a' 

= Y[{l + exp{~(3iE^-f,,))). (A13) 

a 

The occupation number Ua of the single particle state \(pa > for electrons is or 1. It is straightforward to show that 
the statistical operator in Eq. IjAlljl is stationary, i. e. [Ho,po] = and that it fulfills Tr{/5o} = 1. 

3. Expectation values 

Using Eq. PT3)l one finds for the equilibrium mean value (tIq,) of the particle number operator of a single particle 
state \ipa > 

rrS'-U \ X / \ X exp(-/?(£;a - ^,)) . f f \ (ata\ 

irjpoC^CQ'l = daa'{na) = Oaa'- 7 — — "r" = Oaa' JFD(Ea - /is), (A14) 

1 + exp[~l3{Ea - Us)) 

with the Fcrmi-Dirac distribution function 

fFD{E^ - /i.) = — 7-^- ^. (A15) 

l+CXp[l3{Ea - /Us)) 

In agreement with the Landaucr-Buttikcr formalism it then results for the mean value of the particle density operator, 

p(r) = 2Tr{/5o*^(r)*(r)} = 2^/j.z?(i?a-M.)l^a(r-)|' 

a 

/oo 2 
^degs{e)fFD{El+e-fis)\Mr±)fi'^'He,z) , (A16) 

where we establish the continuous hmit by replacing 

E ^ E / '^^■^ ^ E / 9s{e) (A17) 



and the factor two comes from the spin degeneracy. Considering the ansatz Q one can perform the v summation in 
Eq. HA16|I obtaining Eq.©. 



APPENDIX B: HARMONIC PERTURBATION 



In this Appendix we describe our approach to calculate dynamic linear-response properties of open quantum systems 
which are defined in Sect III Al Starting with Eq. H10|l we determine the density matrix for the perturbed system using 
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the von Neumann equation with Hq describing the stationary open system, Eq. (jA10|) . In Hnear approximation one 
finds^ 



Sp{r,uj)= J d^r'no{r,r',uj)SV{r',uj), (Bl) 
where the density-density correlation function (irreducible polarization) is given by 



noir,r',uj) = -J dr cxp{i{uj + ir,)T){[pi{r,T),p{r')])o. (B2) 

The index means the thermodynamic expectation value with respect to the statistical operator pQ as given by Eq. 
(|A11(I and the index / means the operator in the interaction picture. The single particle density operators p are now 
written in the second quantization using the field operators of the scattering states defined in Appendix 1X1 so that 

pi{r,T) = ^ Lp*^{r)ifa'{r) exp [iHot / h]c^^Ca' exp [-iHf)T/h], (B3) 

a, a' 

with a = [sjv) and p{r) ~ pj{r,T = 0). In a standard way one uses the anticommutation relations for the cj^ and Cq, 
to calculate the commutator in Eq. (|B2p and obtains 



no(r,r ,0;) = 2 hm ^ - E^, + + )^«(- )' ^ 

a.cx' 

where we have used Eq. ljA14p for the expectation values. In difference to the usual expression for IIq the Fermi-Dirac 
occupation functions fpD may contain different chemical potentials, either that of the source contact for s = 1 or that 
of the drain contact for s = 2. Furthermore, the single particle energies E^ are given in Eq. IIA10|) and the overall 
factor of two in Eq. HB4|) accounts for the spin degree of freedom. 

In the considered planar structure the potential perturbation depends only on z, 5V{r') = 5V{z). In order to 
develop further Eq. HB4|I . we take the continuous limit [Eq. HA17|l ] and using Eqs. Q and Q one can first integrate 
over r_L obtaining 5{k± — k'±) and after that sum over k±. As a result a z-dependent density modulation Sp{r) ~ Sp{z) 
is obtained which is related to the potential modulation by Eq. Hll|l . and the polarization 110(2;, z', w) is given by Eq. 

®- . . 

Evaluating in the same procedure the perturbation of the expectation value of the z-component of particle current 
density jz(r) = {h/m*)lm I 4'^(r)V2^'(r) ) one obtains 



6jz{\z\<d)= dz'Uo{z,z\uj)SV{z'). (B5) 



The current-density response function for a planar structure has the form flo{z,z',uj) = X^ss' ^f)^ \z,z',uj) with 

<^)(z,z',c.) = lim / de de' ^ , ll^ i-, (B6) 



and 



F^-'\z,z',e,e') = 2-^,.(e)g. (O In | ^ " ^j^ 
27r/5/i^ 1^ 1 + exp [p[ps' - e )J 



h 



2im* 



(V'(^) (e, z)) * ^^(^') (e', z) - V'f^') (e', z) ^ (^(^) (,, ,) 



:(^(^')(e',z'))V'^He,^')- (B7) 



APPENDIX C: LOW FREQUENCY EXPANSION OF THE IRREDUCIBLE POLARIZATION 



For small frequencies we expand in Eq. H12|) for finite rj 
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V V 

FIG. 4: Transformation of the two-dimensional integral l)C4|l . The gray area illustrates the integration domain. 



finding 



and 



/"OO />OC 

pI:'\z,z') = lim / de de 



lim 



de I de 



,F('*")(z,z',e,e') 
e — e' + ihr] 



i v^o Jy^ Jy^ (e - e' + ihr])^ 
We introduce a transformation w = e + e' and u = e' — e so that for a general function /(e, e') 



de / de'f{e,e') = - 
JV, ^ 



/Mo(tl) 
duf(e{v, u), e'{v, u)) , 
-tio(«) 



(C2) 



(C3) 



(C4) 



where we introduce a cut-off energy Voo ^ oo. Furthermore, as illustrated in Fig. 2) uq{v < Vs + Voo) ^ v — 2Vs and 
uo{v > Vs + Voo) = 2V^oo — V. For fixed z, z' we write F^^^^{z, z', e, e') = ax(w, u) with a = 2m* /2tt f3h^ and 



V + u 



N 



N 



V + u 



where 



and 



M{v) = gs{v)(i'^'\v,z)y i;'^'\v,z') 



7V(«) =ln{l + exp {PUis-v))]. 



(C5) 



(C6) 



(C7) 



It is easy to see that xi^i'^) ~ Xii''^^'^) +*X2('1',m) — ^X*{v, —u) so that ~ — and X2(v,ti) 

X2iv, —u). Furthermore, since xi''^^ 0) = one finds for small |?/| the expansion 



X{v,u^ 0) w uxuiv,0) 
where we obtain a real function for the partial derivative with respect to u at u = 



fpD ( 2 ^ ^' 



(C8) 



(C9) 



with the Fermi-Dirac distribution function given by Eq. (HTSJ. This means that in the expansion ljC8|) . the leading 
term in the real part is linear in u while the leading term in imaginary part is parabolic in u. Writing lim,,^o(e — e' + 
ihr])~^ = —PV{l/u) — iTrS{u), where PV denotes the Cauchy principal value and using the symmetry properties of 
the functions xi smd X2, one obtains 



P(^)(z,z') 



dv / du^i^ = / de de' ^ [z,z,e,e) 



2K 



tio('u) 



e- e' 



(CIO) 
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where F*^**^ (z, z', e, e') = F^'^^'^z^z'^e^e') + {z,z',e,€'). In Eq. HC10|I we omitted the principal value operation 

because one obtains in u = a regular integrand due to the expansion in Eq. (|C8|) . Writing^ lim^_^o(£~ + *^^)~^ = 
PV{l/u^) - in{d/du)S{u) it follows from Eq. ^ that 

2 J2V, J-uoiv) " 2 ^214 

-aphn f ^ degUeM^'\e,z)\'\^^-^\e,z')\'fM^- ^is). (CU) 

Here we omitted the principal value operation in the first integral since the integrand is regular at w = because the 
leading order term in X2iv,u 0) is parabolic in u. 
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